Preprint IASSNS-HEP-00/32 
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We re-examine a long-standing problem of a finite-frequency conductivity of a weakly pinned 
two-dimensional classical Wigner crystal. In this system an inhomogeneously broadened absorption 
line (pinning mode) centered at disorder and magnetic field dependent frequency ujp is known to 
appear. We show that the relative linewidth AuJp/cJp of the pinning mode is of the order of one in 
weak magnetic fields, exhibits a power-law decrease in intermediate fields, and eventually saturates 
at a small value in strong magnetic fields. The linewidth narrowing is due to a peculiar mechanism 
of mixing between the stifi'er longitudinal and the softer transverse components of the collective 
excitations. The width of the high-field resonance proves to be related to the density of states in the 
low-frequency tail of the zero-field phonon spectrum. We find a qualitative agreement with recent 
experiments and point out differences from the previous theoretical work on the subject. 



I. INTRODUCTION 

The problem of pinning of an elastic solid by an exter- 
nal random potential has been a long-standing problem 
of condensed matter physics, first addressed by Larkintl 
in 1970. Physical systems where this problem arises in- 
clude charge-densit5j-|Wavesj3 vortex lattices in supercon- 
ductors,H interfacesjj magnetic bubbles,El liquid crystals 
in aerogelsj3 and many others. An interesting example of 
the charge-density wave is a Wigner crystal (WC). Under 
terrestrial conditions a stable WC has been realized in its 
two-dimensional (2D) form,Ll on cryogenic surfaces and in 
semiconductor heterostructures with low carrier density 
and/or in strong magnetic fields .u In this work we exam- 
ine the finite-frequency response of a collectively pinned 
2D WC. The collective pinning is the regime where in- 
dividual defects of the substrate are too weak to signif- 
icantly deform the crystalline lattice, so that the WC is 
well ordered ("unpinned") at small length scales, but the 
cumulative effect of the disorder eventually dominates 
the elasticity at a pinning length, which is much larger 
than the lattice constant. 

The linear-response conductivity aa[3{q,(^) of the WC 
in an external magnetic field is a tensor. We focus 
our study on the real part of its diagonal component 
<7xxi<l — 0,Ci;), the quantity that determines the power 
absorption in the presence of a uniform electric field. If 
the frequency u is not too small, the absorption is domi- 
nated by the collective excitations. In the absence of the 
magnetic field and random potential, a 2D WC possesses 
two branches of such excitations: the longitudinal (L) 
and the transverse (T) phonons. The long-wavelength 
L-phonon is also referred to as plasmon. In a finite mag- 
netic field L- and T-phonons hybridize into the magne- 
tophonons (gapless lower hybrid mode) and the magne- 
toplasmons (gapped higher hybrid mode). It is notewor- 
thy that without the pinning potential, the frequency 



and the oscillator strength of the magnetophonons van- 
ish at g = by virtue of Kohn's theorem. The absorp- 
tion of a spacially uniform ac field is possible only by 
exciting the magnetoplasmon mode, at the cyclotron fre- 
quency ujc = eB/rrieC set by the magnetic field B. The 
pinning shifts the original q — magnetophonon mode 
to a nonzero disorder-dependent frequency Up, broadens 
it, and imparts it with a finite oscillator strength. The 
resulting absorption line, which can be called the pin- 
ning jjiade, was first discussed by Fukuyama and Lee 
(FL).EI'EJ In those early works the linewidth Aujp of the 
pinning mode was predicted to be of the order of iVp itself. 
Strictly speaking, PL considered not a WC but a more 
conventional charge-density wave where the charge mod- 
ulation has only one harmonic, i.e., is cosine-like. Experi- 
ments on such materials, e.g., K0.3M0O3 ("blue bronze"), 
TaSs, etc., performed in zero magnetic field indeed re- 
vealed bipad absorption lines at disorder-depeadent fre- 
quencies.llil A later work of Normand et alxd devoted 
specifically to the WC in a strong B revised some results 
of PL but left unchanged the prediction that Aujp ^>^p- 
At the time, experiments seemed to confirm that.E^I It 
came as a surprise wherLjpiost recent measurements in 
strong magnetic fieldsEdil3 demonstrated that Aujp can 
be more than order of magnitude smaller than ujp. Such 
unexpected findings revived interest in this long-standing 
problem. 

An important step towards resolution of the puzzle 
has been made by Fertig,E3 who showed that long-range 
Coulomb interaction plays an important role in reduc- 
ing the inhomogeneous broadening of the pinning mode. 
However, Pertig used an oversimplified model (commen- 
surate pinning) and did not calculate Aujp directly. His 
results for another quantity may be interpreted as an in- 
dication that Aup/iLip is exponentially small for weak pi, 
ning. More realistic model was studied by Chitra et al 
who considered pinning by a short-range Gaussian ran- 
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dom potential of a general type. We comment on their 
results later in this section. 

In the present paper we study essentially the same 
model as FL and Chitra et al. except we treat the elec- 
trons classically. This limits the applicability of our re- 
sults to the case ^ > Ib, where ^ is th e correl ation length 
of the pinning potential and Ib = y^hc/eB is the mag- 
netic length. The classical approximation enables us to 
focus on the interaction of the WC with disorder, which 
is really the essence of the pinning mode phenomenon. 
We critically re-examine the previous work on the sub- 
ject, identify the physical mechanism responsible for the 
line narrowing, and spell out the conditions under which 
it occurs: 

• the lattice dynamics is inertial (not overdamped) 

• magnetic field is sufficiently strong, such that iVc is 
larger than ujp 

• the pinning is sufficiently weak so that ujp is much 
smaller than the magnetophonon bandwidth 

• the compression modulus A of the WC evaluated at 
the pinning length is larger than its shear modulus 
fi (satisfied for long-range interaction). 

If any of the conditions above is violated, then the line is 
conventionally broad, Aup ~ ujp. 

Although we were unable to find the exact expression 
for Awp, we show that there is an asymptotically exact 
relation between AiOp at intermediate B and the low- 
frequency tail of the phonon spectral function in zero B 
(we remind the reader here that we study classical elec- 
trons, which form the WC at any B). Establishing such 
a connection and elucidating the physical mechanism of 
the line narrowing are our main achievements. 

The phonon spectral function at low frequencies is ex- 
pected to have a power law form 



ImDaaiq = 0, Cj) OC CJ 



2s+l 



(1) 



with a nontrivial exponent s. This additional input and 
further considerations enable us to predict that as lOc 
increases and becomes larger than ujp, Aujp/ujp first de- 
creases in a power-law fashion 

Aujp/ujp (ujp/LOcY (2) 

but then eventually saturates at a value (see Fig. |l|) 

Aup/uop ^ {fx/Xy (3) 

(here A is meant to be ejiaLuated at the pinning length). 
The recent experimentstirtll have probed the high-field 
regime where Eq. (||) is supposed to apply. |_. 

In the aforementioned work of Chitra et aZ.,Efl a set 
of equations was derived, which, when analyzed further, 
also yield Eq. (^), with s = i. This specific value of s 
can be traced down to the fact that Chitra et al. un- 
knowingly rederive a common form of the self-consistent 



Born approximation (SCBA), where Eq. (|l]) with s = ^ 
is satisfied. Since the SCBA is uncontrolled at small w, 
these results are unreliable. Moreover, appealing to the 
known properties of Lifshitz tailsE3 in other disordered 
systems, we argue that the SCBA strongly overestimates 
the low-frequency spectral function. We propose that the 
correct approach should be based on considering certain 
low-propability disorder configurations ( "optimal fluctu- 
ations" ) , which leads to 
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FIG. 1. Relative linewidth of the pinning mode as a 
function of the cyclotron frequency. 

The paper is organized as follows. In Sec. || we discuss 
the ground state of the pinned WC, placing emphasis on 
the parameters that provide the input for the dynam- 
ical problem of interest. In Sec. Ill we formulate the 
model and the general framework for study of the finite- 
frequency response. In Sec. IV we derive Eq. (^. In 
Sec. we study the soft phonon modes at i? = and ob- 
tain the above formula for s. The comparison with the 



experiments is given in Sec. VIl after a brief discussion 
of quantum and finite-temperature effects in Sec. VI. 



II. STATIC PROPERTIES 

The distortions of the lattice in the ground state of a 
collectively pinned Wigner crystal accumulate gradually 
over length scales much larger than the lattice constant 
a. Such distortions can be described in terms of a smooth 
displacement field u'^''^(r). The Hamiltonian of the sys- 
tem H = Hei + Hp is the sum of the elastic term. 



(5) 



H% = ( - ^ i?^(<z) + ^-^Hl{q\ (6) 



q J q 

and the pinning term,0 

Hp = nj C/(r) |^cosK[r-u(°)(r)] - Vu(°) | 
(C/(r)C/(0)) = C{r). 



(7) 

(8) 
(9) 
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The notations here are as follows. The q-integral is 
taken over the Brillouin zone with measure cPq/ {2tt)'^; 
the r-integrals is taken over the area of the system with 
measure cPr. Parameters /i and A represent the shear 
and bulk elastic moduli of the WC, respectively. In 
strong magnetic fields, i.e., at low filling factors they ap- 
proach their classical values J^l /j, « 0.245e^ne^^/K and 
\{q) — —5^ + Y/\q\. Here Ue — 2/\/3a^ is the average 
electron concentration, k is the dielectric constant of the 
medium, and Y = 2'!Te'^n1/ k. The infrared divergence 
of X{q) originates from the long-range Coulomb interac- 
tion. K's are the reciprocal lattice vectors. Finally, the 
correlator C(r) is assumed to rapidly decay at r > ^ 
where ^ < a is the correlation length of the short-range 
Gaussian random potential U{r). 

The static properties of the system defined by Eqs. 
^ can be investigated using classical statistical mechan- 
ics methods. Problems of this type have been studied 
extensively in the past, see Ref. for review and refer- 
ences, and also some recent contributions, Refs. 
p6| . The 2D WC generally conforms to the d — 2, m = 2 
class of elastic media {d is the total number of spacial di- 
mensions and m is the number of components of u) . How- 
ever, some caution is required in transferring the known 
results to the WC case because the L- and T-components 
are not equivalent. The large disparity X:^ fioi the elas- 
tic moduli causes strongly suppression of the longitudinal 
distortions compared to the transverse ones. Therefore, 
in certain formulas one has to use m — 1. In this respect, 
the WC is similar to the Abrikosov vortex lattice where 
the analogous inequality cn 3> cge exists.0 

Let us define the two-point correlation function (rough- 
ness) W{r) = ([u('')(r)-u(°)(0)]2)^/^ Many static prop- 
erties of the system can be deduced from the behavior 
of W{r). Three regimes can be distinguished: random- 
force, random-manifold, and the asympotic (Fig. ph. 







\nr 








random / 


random 


asymptotic 


force / 


manifold 




r / 







FIG. 2. Roughness W{r). 

The random-force or Larkin regime appears for small r 
where W(r) < ^. The last inequality legitimizes expand- 
ing the exponential in Eq. ^ and keeping only the linear 
in terms. The Hamiltonian becomes quadratic, and 
its ground state is easily found to be 



Zi(°)(q) = -zn,5][H0];l(q)X,,^(q + K) 



(10) 



Here and below tildes denote Fourier transforms. 
Averaging over the disorder, we arrive at 

W\r) ^ V'iO) /(I - e^^n [{H"j,)-' + (Hi)-'] (11) 



oc ln(_Rc/^), 



V{z) 



— _ e 

2 



/2 



(12) 
(13) 
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To obtain Eq. (|l|) we cut off the log-divergent integral 
in Eq. ( pT[ ) at q — qc where qc = tt/Rc and Rc is the 
Larkin length defined by the relation W{Rc) — For a 
weak random potential Rc is much larger than a and is 
given by 



(14) 



e estimate for can be obtained in a simpler 
by equating the elastic energy density fiS,^/Rl to 



The 
way,i 

the pinning energy density [C(0)]^/^V^ /R"^. In the last 
expression N = n^R^ is the number of electrons within 
the area Rc x Rc and the square root takes into account 
that individual pinning energies have random signs. 

The elastic-manifold regime is realized at intermedi- 
ate distances, Rc < r < Ra where ^ < W{r) < a. Here 
the roughness grows slower than in the Larkin regime, 
W{r) (X r"^ , < 1. For the 2D Wigner crystal the Flory 
estimate C = I (^s appropriate for the d — 2_m — 2 case) 
should be accurate but perhaps not exact .0 The slower 
growth rate of W{r) is due to a less efficient mechanism 
by which the lattice adjusts to the random potential on 
length scales exceeding Rc. The steepest descent towards 
the nearest local minima of the potential landscape in the 
Larkin regime gets replaced by selection among a set of 
many possible yet roughly equivalent nearby minima. 

The third and the last regime begins when W{r) 
reaches a at the length scale Ra = Rc{a/S,y^''. The ac- 
cumulation of the pinning energy density at such large 
length scales becomes additionally suppressed by lattice 
periodicity effects. Indeed, a uniform shift u —> u -t- a re- 
duces to relabeling the lattice sites with no physical sig- 
nificance. As a result, the roughness grows only logarith- 
mically with distance, W{r) ~ a\n{r / Ra)- The asymp- 
totic In r-roughness corresponds to a "quasi-short-range" 
order: W{r) increases so slowly that it is barely enough 
to remove the algebraic divergencies (Bragg peaks) of the 
static structure factor 



Fiq) 



exp 



i{q-~K)r-^K^W^{r) 



(15) 



K 



at the reciprocal lattice vectors q = K. Remark- 
ably, it also creates enough elastic stress tp-, gener- 
ate unbound dislocations, on the length scalec3 ^jj 

1/2 

_Raexp[Coln ' {Ra/a)]. At this scale u becomes multi- 
valued and Eq. (||) must be modified. 
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Although very interesting, the random-manifold and 
the asymptotic regimes do not play much role in the fur- 
ther discussion. We mentioned them here for the sake of a 
more comprehensive introduction and also to draw atten- 
tion to the difference between two characteristic length 
scales, Rc and Ra- Ra is the length that seems to be 
important only for rather subtle properties mentioned 
above: the destruction of the Bragg peaks and the ap- 
pearance of the topological defects. In contrast, Rc de- 
termines most of the physically important response char- 
acteristics such as the threshold electric field in the dc, 
transport and the frequency ujp of the pinning mode.t2l 
This role of Rc originates from the fact that it is the 
length scale where the dominant contribution ^ fi£^^ /R^ 
to the pinning energy density is accumulated. In this 
sense Rc plays the role of the pinning length mentioned 
in Sec. |. 

Let us now elaborate on the difference between the T- 
and L-correlation functions mentioned above. To do so 
we decompose the elastic displacement into its T- and 
L-components, 



u(q) = qui(q) + [z x q]uT(q), 



(16) 



where q = q/ 
roughness Wa ■ 



q|, and define the two components of 



,(0) 



(r) 



,(0) 



(0)) 



1/2 



a 



L,T. First 



of all, let us examine the short-distance behavior using 
Eq. (IC). For WL{r) we obtain the expression similar to 



Eq. 
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Wl{r)^V'{0) 



o«qr 



27ry2 \a 



(17) 



and so ^^(r) is much smaller than Wrir) — W{r) given 
by Eq. (12). Formula (|l7| ) is certainly valid for r < Rc- 
For larger r we can no longer treat the elastic distor- 
tion |u^°)(r) — u(°^(0)| as small. Nonetheless, we expect 
that the region of validity of Eq. (^^ extends somewhat 
beyond r = Rc. Indeed, instead of doing the perturba- 
tion theory around the ideal lattice state, we can do it 
around the ground state of the incompressible crys- 
tal, A = CXI. In the second approach we have to account 
for a finite pinning energy density, but as long as the 
longitudinal elastic stiffness Yq where q ^ ir/r exceeds 
the characteristic curvature ^q^ of the pinning energy 
density landscape, i.e., for r < Rl — {Y/ ^)Rl ^ R^/a, 
this should be unimportant. Substituting r — Rl into 
Eq. ([l7| ) we find thaLW^i(i?L) <C C> contrary to some 
previous suggestionsEd that Wl{Rl) ~ Wt{Rl) [recall 
that Wt{Rl) ^ W{Rl) » W{Rc) = 

The large-distance (r > Rl) behavior of Whir) is un- 
accessible within the perturbation theory. However, we 
should get qualitatively correct results by .using the Gaus- 
sian variational replica method (GVRM).Eij Assuming for 
simplicity that Rl ^ Ra and suitably modifying GVRM 
for the present problem, we find 



Wl{r) oc V^'(O) 



1 — cos(qr) 



HO{q)[Hliq)+na 



(18) 



where Ila = f^Ra^- It is easy to see that WL{r) tends 
to a finite limiting value as r ^ oo. So, the long-range 
positional order is absent in the T- but is-^reserved in 
the L-part of the elastic displacement field.EZi 

To conclude the review of the statics, let us consider a 
matrix of the second derivatives of Hp, 



S(r) = 



au(o)au(o) 



= 7ie^KK;7(r) cosK[r- u(°)(r)]. 

(19) 



K 



Matrix S plays an important role in the dynamics and 
we need to understand its properties. 

In the case of a short-range disorder we study here, 
the correlations in U(r) decay more rapidly with r than 
those of u(*'^(r). This enables us to find the two-point 
correlator of Safd with little effort. It suffices to use the 
approximation u(°)(r) = in Eq. (|l9|), which leads to 



(5Q^(ki)S'^^(k2)) - {Sap.){S„f3) 



V"iO) 



X (<5q^<5i//3 + SapS^„ + (5Qi/5^/3)(5kik2^^ 



(20) 



for the WC with the hexagonal lattice and in the long- 
wavelength limit fci, fc2 *C l/o- Here is the area of the 
system. Another important property of Sap is a nonzero 
mean, 



(21) 



Parameter is positive, in accordance with the pinning 
phenomenon: the crystal in its ground state is distorted 
in such a way that the electrons are preferentially situ- 
ated near the minima of the random potential U(r) where 
its curvature is positive. Sq can be estimated to be 



So 



V"{0) /Rc 
27r/i 



(22) 



using the same procedure as for deriving Eq. ([T^). As- 
suming that |u*-°-'(r) — u(°^(0)| is "small," we expand the 
exponential in Eq. ( [l^ ) to obtain 

5„^(q) ^-'-^Y. KcKpK^ I <°)(qi)C/(q - qi - K). 

(23) 

Combining Eqs. (^^ and (|2^), averaging over the disor- 
der, and using the inequality H^{q) ^ H^{q), we get 



5o -V"{0) / (q)] 



(24) 
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which leads to Eq. (|2^). Here we again had to cut off 
the infrared logarithmic divergence by hand, a,t q = qc = 
tt/Rc- More sophisticated approximation schemes such 
as the GVRM would implement this cutoff more grace- 
fully but as long as we are not interested in the numerical 
factor in the argument of the log, the result is the same. 

At this point we conclude the discussion of the ground 
state properties of the pinned WC, as we are now ready 
to address its dynamical response. 



III. FINITE-FREQUENCY RESPONSE: 
PERTURBATION THEORY AND QUALITATIVE 
CONSIDERATIONS 

To study the dynamics of the WC we introduce a 
time-dependent displacement field u, which is the de- 
viation of the total displacement from its ground-state 
value u'^°) (r) . We will restrict ourselves to the harmonic 
approximation where the action is quadratic in u. The 
response of a harmonic oscillator is the same in quan- 
tum and classical mechanics, so we can use the conve- 
nient imaginary-time formalism without compromising 
our original intention to treat the system classically. The 
resultant action contains two terms: one which describes 
a uniformly pinned WC, 



An 



1 1 

2h^ 



I ut(q,za;„)[D°]-iu(q, 



iuj„), (25) 



and the other which takes into account fluctuations in 
local pinning strength, 

qi q2 

(26) 

Here /? — \/{kBT) is the inverse temperature, w„ — 
2Tin/{h(3) are the bosonic Matsubara frequencies, = 
S — 5oI, I is the identity matrix, D" = D°(q, iw„) 
is the phonon propagator of a uniformly pinned WC 
(cf. Refs. ^Ol and p2|) 



f ^0 + 



Hi + So 



(27) 



p = rriene is the average mass density, and Rq is the 
0(2) rotation by angle a,rg{qx + iqy)- We also define the 
disorder-averaged propagator D(q, iun), 



D(q, iujn) 



J I?ul?u*u(q, z^„)u*(q, iu!n)e~ 



(28) 



where A = Aq+Ai and is again the area of the system. 



The quantity we set out to calculate is the ae conduc- 
tivity cr{q,uj), which in this model is given bjnlj 



(T{q,Lj) = — ie^rigCj D(g,ia;„) 



(29) 



The conductivity can also be expressed in terms of the 
phonon self-energy 



n(q,c^) = 5oI+ (D-i-[D"]- 



(30) 



We are interested primarily in the case q = where the 
most general form of n(q, w) consistent with rotational 
symmetry is 



UapiO, Uj) ^ (5a/3H(w) - iSappLULOcfxyii^), 



(31) 



fxy{^) being the relative correction /S.pxy/ p'^y to the bare 
Hall resistivity p^yiuj) = B/{neec). Approximate calcu- 
lations presented in this paper give fxy — 0, which leads 
us to believe that fxy{Lij) must be small for weak pinning. 
We choose to-aeglect such fine details and to assume that 
fxy vanishesHS In this case 11(0, uj) is a scalar and (Txx{'^) 
is given by 



Re <7xx {'^) — e^TigCj Im ■ 



H(0,e) 



[n(0,e)-e]2~«,' 
where we introduced convenient "energy" variables 



puj 



(32) 



(33) 



From Eq. (|3^) one can see that in strong magnetic fields, 
LUc 3> iOp, the power absorption in the uniform electric 
field takes place mainly within the frequency interval of 
width 



Alu.„ 



ImH(0,ep) 



puc 



centered at w = Wp, where 



Re H(0, ep 

pUJc 



(34) 



(35) 



The last equation is the implicit definition of ep. In the 
following we will work predominantly with the "energy" 
variable e rather than the frequency uj. 

At small e the calculation of 11 is a difficult problem 
but at large energies, e 3> Cp, the first Born approxima- 
tion suffices. The only parameters needed to implement 
it are the mean and the variance of the matrix elements 
of S given by Eqs. @) and (|l|). For g < 1/a we obtain 
that HQ^(q, e) = 5a^H(e), where 



H(e) - 5o + r'(0) J trD"(k,e). 



(36) 



H has both imaginary and real parts, corresponding to 
the broadening and the frequency shift of the pinning 
mode, Eqs. (|34) and ([3^), respectively. The imaginary 
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part comes solely from the pole(s) of D''(k, e), i.e., the 
solutions of [iJ^(k) + So - e][^°(k) + 5o - e] - ee^ = 0. 
A more accurate expression for 11 can be obtained using 
the self- consistent Born approximation (SCBA), 

n(e) = 5o + V"{0) J trD(k,e), (37) 



as long as the full propagator D(k, e) has a pole at 
|fc| S> Qc — T^/Rc- Under this condition the diagrams with 
intersecting lines not included into the SCBA are sup- 
pressed, much like they are suppressed in a dirty metal 
with kp ^ kp being the Fermi momentum and I 
being the mean free path. The indicated condition is 
satisfied when 



e 3> min < Hi 



^3 



(38) 



Here and below A is meant to be evaluated at q — qc un- 
less it is indicated otherwise. Parameter Ho, which can 
be estimated to be 



Ho ^J.q^ 



(39) 



represents the real part of the self-energy at the lowest 
e allowed by inequality (|38|). Ho is smaller than So by 
the logarithmic factor due to the partial cancellation be- 
tween the first and the second terms in Eq. (37). This 
unfortunate cancellation leaves us only with the order of 
magnitude estimate (p9|). 




(A/yu)n„ 



FIG. 3. Top: Pinning mode location as a function of 
the magnetic field. Bottom: SCBA domain of applicabil- 
ity. 

When inequality ( |3^ ) is violated, the SCBA becomes 
an uncontrolled approximation. The only SCBA result, 
which can presumably be trusted at such e, is a slow de- 
pendence of Ren(e) on en ergy . This is because the real 
part of the integral in Eq. (B^) is dominated by large fc's 



where the lowest-order perturbation theory is valid. ..Tp 
estimate Cp we can take Ren(e) « Hq, which yieldsE2lE3 



Cp ^ mm <^ Ho, — 



(40) 



represented graphically in Fig. y together with the SCBA 
domain. 

As one can see from this graph, for Cc < Hq (weak 
magnetic fields) there is no parametric separation be- 
tween the SCBA domain and the curve e = ep] therefore, 
the SCBA is qualitatively correct even along this curve 



-Iln. Hence, Aujr, 



UJr, 



where it predicts Imn(ep) 
The resultant broad absorption maximum is sketched in 
Fig. |. 



Recr, 




FIG. 4. Frequency dependence of the conductivity in 
weak magnetic fields. 

If A and fi are comparable, the same argument works 
also in strong magnetic fields (cc » Ho) where axx{i^) 
looks qualitatively similar, except the position of the 



maximum depends on the magnetic field: ujp = LOpQ/tUc, 
where Wpo is the pinning frequency at i? = 0. On the 
other hand, if A ^ /x, as for the WC, the pinning mode 
is situated in the interior of a parametrically wide range 
of £ where the SCBA fails. Thus, the calculation of Aujp 
requires other methods. Later we will show that Aujp in 
strong magnetic fields is related to the low-frequency tail 
oi CTxxi'-^) in the absence of the magnetic field. Although 
calculating the functional form of such a tail is another 
nontrivial problem (after all, it is beyond the SCBA!), 
this intriguing relation is sufficient to establish that the 
absorption-line, narrows down, in agreement with the ex- 
periments .tjil3 

We unfold our argument gradually over the remaining 
sections. In this section we start implementing this task 
by clarifying why Uxxii^) is virtually independent of the 
longitudinal stiffness A provided A ^ /i and 5 = 0. 

Since action A is quadratic in u^, the L-phonon de- 
grees of freedom can be easily integrated out, leading to 
the effective Hamiltonian for the T-phonons, 



H7 



(H 



eI)-iS 



(41) 



where H^, H^,, and should be understood as oper- 
ators. The first two are diagonal in the basis of plane 
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waves, with matrix elements H^{q) and H^{q), respec- 
tively; Si's have both diagonal and off-diagonal matrix 
elements: 



'S'rqq' = (^a;3trS — Sa[3)qaq'i3, 



(42a) 
(42b) 
(42c) 



Finally, Sx+ = Sx + ^/ee^I- Integrating out is an 
exact algebraic transformation, which preserves the spec- 
trum of the collective modes: the resolvent Dxi^) of the 
operator Ht — el has poles at the same e as the full prop- 
agator. In fact, Dt, which can be written in the form 



Dt = ([DJ]"i - S* 



where 



(HO 



el) 



S* = St,,(H^ + Si-eI)-iSx+, 



(43) 

(44) 
(45) 



is nothing else than the T-T component of the full 
phonon propagator before the disorder averaging. The 
self-energy of the averaged propagator at g = is a scalar 
(see above); therefore, the T-T component determines 
the entire propagator for this particular q. 

Let us now show that in zero magnetic field the mixing 
between the T- and L-phonons represented by S* has lit- 
tle effect on D^. Indeed, in the diagrammatic expansion 
of (Dr(e)) in powers of SSt and S*, the typical momen- 
tum transfer k for e ^ IIq is of the order of gc- A single 
occurence of S* contributes k^\V"{0)\/\P ^ i^J./X)UQ 
to the self-energy, compared to ~ Hq from SSt- Dia- 
grams with multiple occurences of S* are suppressed by 
even higher powers of the small parameter /i/A. Thus, in 
the parameter range relevant for the observation of the 
pinning mode 



(DT)(q,e)^ (DJ)(q,e)^ 



1 



i/^(g)+n*(q, e)-e' 



(46) 



Note that is determined only by the shear modulus 
fi and the disorder in the T-T channel SSt, while bulk 
modulus A drops out. It follows from this discussion that 
for A ^ /z we can calculate the g = response pretending 
that the L-degree of freedom does not exist. 

Let us indeed imagine that the L-phonons are forbid- 
den and try to investigate the nature of the T-phonon 
eigenstates that would compose the pinning mode, i.e., 
the states that would respond to the uniform electric field 
(still at S = 0) . To do so we need to analyze the solutions 
of the eigenvalue problem for the single-particle Hamil- 
tonian H* = H^. + S^. In the uniformly pinned WC, 
St = ^qI and the eigenstates are just the plane waves 
labelled by momenta k of the Brillouin zone. In the ac- 
tual random system the eigenstates are wavepackets of 
the plane waves with characteristic spread of momenta of 
the order of L^^, where Lt ^ wt/ |IinII| is the T-phonon 



mean free path, vt = ^e/I playing the role T-phonon "ve- 
locity." Since |Imn(ep)| ~ t 



ml, 



at e ~ Cr, we have 



Lt ^ Rc- Hence, the eigenstates that respond apprecia- 
bly to a (7 = external drive are wavepackets built from 
< k < qc plane waves. For such wavepackets the aver- 
age k is of the order o£-the inverse mean free path and 
the loffe-Regel criterionc3 maintains that these states are 
localized. In other words, Rc is not only the mean free 
path but also the localization length of the states within 
the pinning mode. 

Let us define the B = T-phonon density of states. 



Ke) 




1 



;DT)(k,e). (47) 



At e 3> ep it tends to the bare value vq = (47r/i)^^, and at 
e ~ ep it is only slightly smaller. It is easy to see then that 
within an area Rc x R^, and energy interval < e < Cp, 
there is typically only one (localized) state. The picture 
that emerges is illustrated in Fig. |^: we have a collection 
of localized states of roughly the same size (i.e., the in- 
verse participation ratio) and roughly the same distance 
from each other in real space, both of the order of Rc- 
The broad distribution of shapes and sizes characteris- 
tic of localized states in disordered systems yields a large 
inhomogeneous broadening Awp ~ cOp of the absorption 
line, as we found earlier based on a different line of rea- 
soning. 

In early works on pinning a ieuristic imagery of "do- 
mains" was sometimes invokedEHl to describe the ground 
state structure. However, it was never clear where to 
draw the boundaries between different domains. This 
difficulty is partially resolved in our picture of phonon 
localization where domains can be defined as areas where 
individual localized phonon wavefunctions are appreacia- 
ble. However, modern understanding of the subject re- 
viewed in the previous section no longer appeals to any 
"domains," and so this interpretation may be regarded 
as a historic sentiment. 




FIG. 5. Sketcli of spacial distribution of T-plionon 
states contributing to the pinning mode at _B = 0. 

One, more historic comment is in order. Previous 
worked on phonon localization found that the phonon 
eigenstates remain extended even in the w — > limit, 
seemingly in contradiction to our results. In fact, there 
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is no contradiction because the model we study here and 
the conventional formulations, where disorder originates 
from the defects of the crystalline lattice, e.g., substitu- 
tions by impurity atoms, are quite different. While in 
those conventional models the defects move with the lat- 
tice and their influence rapidly diminishes at small k's, 
the pinning potential for the WC is static, so it is not 
affected by the lattice motion, ft couples to the phonons 
via the matrix S whose fluctuations have approximately 
the same rms magnitude at all k's within the Brillouin 
zone, see Eq. (|2^). This is the reason why in our case 
the disorder effects are much stronger, causing the local- 
ization of the phonon eigenstates. Actually, our model 
has more in common with those of noninteracting elec- 
trons in dirty metals and semiconductors. (At the one- 
particle level the distinction between the bosonic statis- 
tics of phonons and the fermionic statistics of electrons 
is irrelevant). 

Concluding this section, we would like to reiterate one 
of its main results, that the presence (or absence) of the 
stiff L-phonon degree of freedom affects the B = 0, q = 
response only weakly. Naturally, it makes the fate of the 
L-phonons seem quite mysterious. The finite-_B dynam- 
ical response has also been barely touched upon. These 
gaps in understanding will be filled in the next section. 



IV. SCATTERING OF L-PHONONS AND 
PINNING MODE LINEWIDTH IN A FINITE 
MAGNETIC FIELD 

In the previous section we discussed the scattering and 
localization of the T-component of WC lattice vibrations. 
Let us now turn to the L-component. This time we inte- 
grate out ut to obtain the effective Hamiltonian for the 
L-phonons, the corresponding Green's function, and the 
self-energy: 



Hi = + Sl 
T>L = (Hi - el) 



(H^ 



el)-isl 



el). 



(48) 
(49) 
(50) 



Since the self-energy n(q, e) of the original system (be- 
fore integrating out ut) is a scalar at q = 0, 11^ must 
obey the relation 



ni(0,e) =n(0,e) 



n(0,e)-e- 



(51) 



Our first task is to reproduce the results of Sec. Q by 
demonstrating that in sufficiently weak magnetic fields 
IIl indeed has the above form with 11 ~ 11*. 

There are two types of processes that contribute to 11^: 
the intraband scattering (L-L channel) and the scattering 
through an intermediate T-state (L-T-L channel). They 
originate from the second and the third terms of H^, re- 
spectively. The intraband scattering will be treated in 
more detail in Appendix ^ where using a combination 



of the renormalization group and instanton methods we 
show that it is exponentially small. For our current pur- 
poses a weaker result is sufficient, that the intraband con- 
tribution IIl-^l to Hl is much smaller than Ho. This can 
be established without detailed calculations. Indeed, the 
effective random potential in the L-L channel has ba- 
sically the same mean and rms fluctiations as St; there- 
fore, Sl and St can be regarded as disorder of approx- 
imately the same strength. Then the alluded inequality 
IIIl^lI <C |n| ~ Ho follows simply from the fact that 
L-phonons have much steeper bare dispersion relation 
e{q) = X{q)q'^ than T-phonons and are scattered much 
weaker by disorder of the same strength. It is legitimate 
to neglect the L-L scattering altogether by replacing S l 
with its average value ^qI, and to concentrate exclusively 
on the L-T-L processes. In this case, much likepfor a 
dirty semiconductors with two bands of carriers ,lj the 
diagrammatic expansion of 11^ can be formulated as the 
sum of all one-particle irreducible (IPI) graphs involving 
Sx+ and Dt, see Eqs. (&E5|) and Fig. |a-b. 



(a) 



(b) 



(c) 



(d) 



IPI <e 



Dt 



n, = - 



d; , d; 



-o o- 



=® — o 



FIG. 6. (a) Diagrams for II*. Filled circles and solid 
lines represent SSt and D^ = (H^' -I- Sol — el)~^, respec- 
tively. Dashed lines symbolize contractions, (b) Diagram- 
matic definitions of Dt and Hl- (c) Diagrams without S*. 
Diamonds represent factors i^/ee^. (d) Some representive 
diagrams for Hl not included in (c). 

The crucial point substantiated below is that S* can be 
neglected, i.e., that Dt can be replaced by D^-, provided 
A S> /i and the magnetic field is not too strong, 

ee < {X/fi)Uo. (52) 

For such Ec only two graphs contribute to IIt (Fig. ^), 
which is a significant simplification. However, we still 
have to specify how to average over the disorder. Naive 
averaging of Sx and separately from each other yields 

nT(q, e)^So~ eec(D?,)(q, e) + V"{0) J (D^)(k, e) 

k 

(53) 

but it is certainly not accurate because Sx and i5St are 
not independent. Indeed, their two-point correlator 



8 



(S'xqiqS'Tqsqs) « (qiq2)(q3q) + (qiq3)(q2q) 

+ (qiq)(q3q2) (54) 

is in general nonvanishing. 

In the special case of q = a further progress is possi- 
ble. Comparing the above equation with 

(S'Tqiq^'Tqaqa) « (qiq2)[q X qsjz + (qiq3)[q X q2]z 

+ (qiq)[q3 X q2]z, (55) 

we observe that the former is transformed into the latter 
by the replacement q ^ [z x q] . If the incoming momen- 
tum q is zero, rotating it by 7r/2 has no effect. Thus, for 
q = we can replace Sx in the second diagram of Fig. ||c 
by SSt, after which it becomes identical to the totality 
of diagrams in Fig. |a for n*; therefore, 



ni(0,e) = n*(0,e) 



n*(o,e) 



(56) 



As explained in Sec. pll| , 11* ~ 11; hence, we will succeed 
in reproducing Eq. (pl|) as soon as we show that inequal- 
ity ( |52| ) is indeed the relevant condition for dropping S* 
in the diagrammatic series for 11^. To do so we develop 
one step further the argument of Sec. III. One-particle 
irreducible diagrams containing S* have at least one L- 
phonon line with typical momentum transfer k ^ Qc, see 
Fig. ^d. If Ec = such diagrams are suppressed by at 



least a factor of fi/X, as we found previously in Sec. [I] 



However, if €c ^ 0, these diagrams also generate terms of 
the order of 



1 



Ml 



(57) 



which act as an additional self-energy correction to 
in Fig. They can be approximately accounted for 
if we replace e in the argument of 11* in Eq. (56) by a 
renormalized value e, 



ni(0,e) =n*(o,e~) 



n*(o,?) 



Mc 



(58) 
(59) 



The difference between e and e can be safely ignored 
as long as the inequality ( p2[ ) is satisfied, in which 
case Eq. ( |5^ ) is asymptotically exact. In stronger fields 
Eq. (|5^) applies, which gives only the order of magnitude 
of Hi because of the uncertainty in the parameter Ci. 

The meaning of Eqs. ( |5^ ) and (58) can be elucidated 
returning to the the original formulation where both T- 
and L-degrees of freedom are present: 

n(0,e) = n*(0,e), ec«(A//i)no, (60a) 
= n*(0,e), ec»(A/M)no. (60b) 

We see that 11(0, e) remains magnetic field independent 
up to rather high B. As explained above, this occurs 
because the enhancement of the L-T mixing due to the 



Lorentz force is strongly impeded by the disparity of the 
two elastic moduli. 

Let us now evaluate the consequences of the obtained 
equations (30a) and ( |60b| ). In the intermediate field 
range, IIq ^ Ec ^ (A/^)no, where Eq. (60a) applies, 
the consequences are two-fold: the suppression of ep, 
which is well known, and the suppression of the pinning 
mode linewidth, which is nontrivial. Indeed, since 
becomes much less than IIq [Eq. (^0[)], it slips into the 
low-energy tail of the zero-field phonon spectrum where 
Imll ^ IIq; hence, Acijp ^ ujp. The connection between 
the linewidth of the absorption line in strong magnetic 
fields and the low-frequency phonon modes in zero mag- 
netic field, which we just established, is the keystone of 
the present paper. 

The properties of the zero-field soft modes will be dis- 
cussed in more detail in Sec. [v| where we argue that they 
give rise to the power-law dependence of Imll on e. 



Im^(e)~-C2^o(e/^o)^ e«no. 



(61) 



with exponent 5 = 3/2 and numerical prefactor C2 ^ 1. 
Combining Eqs. (||-§|) and (^^-§l|), we find for the 
intermediate- i? regime: 



/l(^) 



l-«/i(^) 



const, UpQ <C a; <C Wc, 



LUpO < tJc < LOpoy/X/fJ,. 



(62) 

(63) 
(64) 



In this regime the absorption line narrows down accord- 
ing to Eq. (|^), which we reproduce here for convenience: 



(65) 



The high- field regime is described by Eq. (30b). In this 
case the relative linewidth is field-independent and is of 
the order of 



(66) 



as given by Eq. (||) and illustrated in Fig. |l|. Al- 
thou gh E q. ( |60b[ ) was derived with much less rigor than 
Eq. (60a), the saturation of Aujp/ujp in strong fields 
is certainly to be expected. Indeed, in strong fields 
the dynamics is dominated by the Lorentz force; there- 
fore, e must enter in the combination eec, not by itself. 
But then a phonon eigenstate, which for a given (large) 
ec = e* = p{^*)'^ has an energy e^, is also an eigenstate 
of the system for larger Ec, with the eigenvalue ei(e*/ec). 
Hence, as the magnetic field increases, all relevant eigen- 
frequencies scale inversely proportional to Uc, while the 
conductivity varies in the self-similar way. 



(67) 
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and posseses constant Aujp/uj^ The explicit expression 
for crxxii^) is similar to Eq. (pSi), 



1 - if2{Lu) 



const, r2 ^ ^ LJc, 



= C3 a/A///, Wc>WpoVA//i, 



(68) 
(69) 
(70) 



where C3 ~ 1. 

Formulas ( |63| ) and ( |69| ) are our final results. They de- 
scribe the entire lineshape of the pinning mode both in 
intermediate and strong magnetic fields. However, their 
derivation was presented in diagrammatic rather than 
physical terms. Next, we will give an alternative deriva- 
tion, which elucidates the physics of the line narrowing 
and also helps to clarify the structure of the phonon 
eigenstates that compose the pinning mode. 

We will start with the qualitative picture of localized 
T-phonons developed in Sec. Ill and add a new ingre- 
dient, the stiff L-degree of freedom. It is clear that an 
admixture of the L-component produced by a joint ac- 
tion of the disorder and the Lorentz force makes phonon 
eigenstates much more extended in real space. Unlike the 
softer T-phonons, the stiffer L-ones cannot be confined in 
small areas of size Re- In order to understand the large 
scale structure of the eigenstates, we can coarse-grain 
the system by integrating out the degrees of freedom on 
the spacial scales between a and Re- An insight on the 
form of the effective Hamiltonian H'j^^ after the coarse- 
graining is given by the spectral decomposition of the 
matrix element 



1 



HO 



St — el 



r = 



(71) 



Here and UTi{v) are the eigenvalues and the eigen- 
functions of localized T-phonons. Since each of uxi is 
localized within an area of size i?c, we expect that after 
the coarse-graining the numerators UTi{'r)u^^{r') trans- 
form into local operators R~^S{r — r')S{r — r^), where 
is the center-of-gravity of the ith mode. As for the de- 
nominators, the discussion in Sec. [II shows that in weak 



and intermediate magnetic fields [defined by the inequal- 
ity (|5^], Ei's should remain unaffected. Finally, naive 
averaging of in Eq. (^ yields S'qI but the interaction 
with high-fc T-modes, which can be calculated within the 
SCBA-type perturbation theory, renormalizes Sq to Eg. 
The resultant coarse-grained Hamiltonian is 



H 



etc 



RlSiv~n), (72) 



where ^ 1 are random and € (0, 2no). Without 
losing essential physics, we can assume that r^'s form a 
regular square lattice with lattice constant Re- 



What are the properties of the obtained lattice model? 
One of them is the Lorentzian tails of the distribu- 
tion function P{U) of the on-site disorder terms Ui = 
(c^Hq -I- eec)/(e — e^). For example, if e ^ fp, then 



PiU) 



JJ2 ' 



\u\ 



(73) 



where v{£} is the T-phonon density of states introduced 



in Sec. Ill 



Clearly, the variance (C/^) of the on-site 
disorder is unbounded. One can therefore expect much 
stronger disorder effects than in the case of a Gaussian 
distribution with the same typical values ofUi, i.e., same 
(In |C/|). Let Ll be the mean- free path of L-component of 
the low-energ5{-^honons. As explained in Sec. HI, loffe- 
Regel criterioncJ suggests that Ll is simultaneously their 
localization length. For weak Gaussian disorder and un- 
screened Coulomb interactions is exponentially large, 
see Appendix Now we will show that in the presence 
of the long tails ( |73| ) , is only as a power-law function 
of the disorder strength. The crucial point is that the 
scattering is dominated by a few largest C/^'s. As soon 
as we realize this, we can replace P{U) by the Cauchy 



distribution 



PiU) 



1 



TT ([/ - C/o)2 + r2 ' 



r(e) ^ n'^RMe) (74) 



on the grounds that it has the same tails. At this point we 
make a reasonable assumption that correlations among 
Ui^s at different sites can be neglected, and arrive at the 
famous LloycLmodel. The self-energy can now be calcu- 
lated exactly,E3 



HL(e) =no-{/o(e)-zr(e). 



(75) 



while the mean-free path is the solution of the equation 
ijO(i-i) ^ jlmHil, which gives Ll ~ Y/T. 

To compare with our earlier results we need to know 
V and Uo- The integral in Eq. ( ^7[ ) is determined by 
|k| < qc, therefore, 



^(g) TT2 ^ ^ 

^^0 



Hr 



e<no. (76) 



The estimate of Uq is Uq ^ eec/Ilo- Combining these 
together, we find that Eq. ( [75| ) is in agreement with 
Eq. (|5^), which validates our mapping onto the Lloyd 
model. 

To describe the dominant scattering mechanism in the 
Lloyd model notice that Im H^ comes from the inhomo- 
geneous broadening of the phonon energies, i.e., the en- 
ergy shifts 5e caused by on-site disorder Ui. For each 
localized L-phonon the largest energy shift comes from a 
site with largest Ui within the localization area Ll x Ll, 
which is typically Ui ~ (Ll/Rc)^^ for the Cauchy distri- 
bution (74). Since the amplitude of the L-phonon wave- 
function at the position of the oscillator is ml '--^ 1/Ll, 
the resultant energy shift is <5e ~ Uiu\R^ ~ F. The 
energy shifts from other sites are subdominant; thus, 
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ImllL r, in agreement with the exact result, Eq. (75). 
In the original formulation large Ui come from the local- 
ized T-oscillators whose energy is almost in resonance 
with e, i.e., je^ — e| ^ l/i^{e)Lj^. A few of such resonant 
scatterers will localize the L-component on the scale of 
the mean- free path L^. The proposed picture resembles 
the situation in a crystal pinned by strong dilute impu- 
rities, the role of impurities played by the resonant sites. 

The structure of a typical phonon is illustrated in 
Fig. R It involves of the order of {Ll/Rc)'^ T-modes of 
Fig. pFoscillating with the same frequency imposed by the 
mutual coupling mediated by the L-phonon. The oscilla- 
tion amplitude is roughly the same for all the T-modes 
except for a few resonant ones, where it is much larger. 
These resonant sites "drain" the energy of the L-degree 
of freedom thereby localizing it. 




OOQOOOOOOOOQO 



FIG. 7. Schematic structure of the localized phonon 
eigenstate in the crystal with both T- and L-degrees of 
freedom present. Circles represent the T-modes of Fig. ^, 
the resonant modes are shaded. 



Let US now describe the evolution of the phonon eigen- 
states as the magnetic field increases. Since the pin- 
ning mode frequency goes down, lower and lower fre- 
quency resonant sites are required to effectively scatter 
and localize the mixed L-T phonon vibrations. Such sites 
are therefore the soft modes, the oscillators of unusually 
low frequency, which shape up the tails of the zero-field 
phonon spectrum (Lifshitz tails). The soft modes ap- 
pear due to very rare, peculiar disorder configurations. 
The lower the frequency, the mode dilute in real space 
they are. Correspondingly, the mean-free path (local- 
ization length) of the magnetophonons becomes larger 
and larger. The narrowing of the absorption line is then 
similar to the motional narrowing phenomenon. Even- 
tually, in very strong magnetic fields, the enhancement 
of the L-T mixing by the Lorentz force start to dimin- 
ish the frequencies of the localized T-phonons. At this 
point slipping of the pinning frequency Up deeper into the 
soft-mode tail stops. The wavefunctions of the phonons 
cease to change, only their eigenfrequencies continue de- 
creasing in inverse proportion to the magnetic field. As 
a result, Auip/tUp remains constant. 



The physical structure of the pinned phonon modes 
being clarified, the only important questions that remain 
to be considered are the nature of the soft modes and 
the derivation of Eq. (|6l|). This will be the subject of the 
next section. 



V. SOFT MODES 



It is clear that the soft phonon modes must come from 
rare places where pinning is unusually weak. The actual 
calculation of the density of states ^{e) in the low-e tail 
is however nontrivial. The method most suitable for the 
task seems to be the method of optimal fluctuation. It 
has been successfully employed for calcttlation of the Lif- 
shitz tails in other disordered systemsJHJ where it proves 
to be asymptotically exact. The method is based on the 
following idea. The energy e of a given eigenstate is de- 
termined in a complicated way by the distribution of the 
random potential within the entire area supporting the 
eigenstate and so there are many different random po- 
tential configurations, which give the same eigenenergy. 
If e is small, any such configuration is untypical, i.e., a 
fluctuation of some sort. The method is based on the 
assumption that certain fluctuations have much higher 
probability than the rest and dominate the quantity of 
interest, i^(e) in this case. The objective is then to design 
such an optimal fluctuation and evaluate its probability. 

Let us mention one type of fluctuation, which is not op- 
timal. This is a configuration where the random potential 
is stro ngly suppressed in a large area L x L, where L ^ 
■\//i/e. (It may be useful to recall at this point that we are 
investigating the soft modes of a system where L-phonons 
are integrated out) . The probability of such a fluctuation 
can be estimated by mentally dividing the area into un- 
correlated blocks of size Rc x Rc and multiplying together 
the probabilities ^ e/IIo that the random potential is 
suppressed within each block. The total probability is 
exponentially small ~ exp[— (Ilo/e) ln(no/e)]. This is a 
general situation for fluctuations spread over a large area 
in real space. Thus, the optimal fluctuation must be of 
the smallest possible size. It is easy to see that this size 
is Rc- Indeed, to get a small phonon eigenenergy e in a 
small volume, the positive "kinetic energy" ^L~^ must 
be accurately compensated by the negative "potential en- 
ergy" {5St)l- The fluctuations of the latter averaged 
over the area L x L are of the order of Ilo{Rc/ L). Thus, 
for L ^ Rc the cancellation of fj,L~^ can be done by a 
typical fluctuation of SSt and the probability of such an 
event would not have any exponential suppression. In- 
stead, we expect the power-law dependence of i^(e) on e, 
Eq. (g). 

In order to possess a soft mode, the system must have a 
soft direction, i.e., a certain collective coordinate X such 
that the total energy of the system E as a. function of X 
has a very shallow minimum. Let A" = be the ground 
state. In its vicinity E should be Taylor-expandable, 
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E{X) = ^X'^ + pX^ + -fX* + ... (77) 

The coefficient a has the meaning of an effective spring 
constant of the local oscillator, while the phonon energy 
e is determined by the ratio of the spring constant and 
the effective mass. The latter is proportional to the area 
involved in the oscillations. As we argued above, this 
area is of the order of for all e < e^. Therefore, a 
should scale linearly with e in the limit e — > 0, 

a = de. (78) 

If we find a way to estimate the probability of such untyp- 
ically small a, we will succeed in calculating the exponent 
s of the power-law function 

This kind of calculation is hardly possible without 
specifying X. Actually, ascribing the precise meaning to 
X amounts more or less to designing the optimal fluc- 
tuation. At the moment we do not have a good un- 
derstanding how to do it in the 2D case. However, in 
one dimension the general structure of the calculation is 
clear. We are guided by the earlier work on thft-aubject 
by Feigelman et alSB and by Aleiner and RuzinEj 

In the case of a weakly pinned ID elastic chain, the 
role of the desired collective coordinate X can be played 
by the elastic displacement field u at an arbitrary point 
in the middle of the chain. This point divides the sys- 
tem into two parts, which communicate only through the 
single variable u. For any fixed u serving as a boundary 
condition, we can find separately the ground state ener- 
gies i?< (u) and £'> (u) of the left and the right halves of 
the chain. The ground state of the whole system is de- 
termined by minimizing the sum E{u) = E^{u) + Ey{u) 
with respect to u. Both E^{u) and Ey{u) are periodic 
function of u with the period equal to the lattice constant 
a. On their period they have several (typically, of the or- 
der of a/^) minim|ar,and maxima. The beautiful idea of 
Aleiner and RuzinEJ (totally overlooked in Ref. ^3|) was 
that the most efficient way to generate a soft mode is via 
a frustration, when a maximum of E^ (u) occurs near a 
minimun of Ey{u). In this case the second derivative of 
E, i.e., a — E" , can be very small even though each of 
the second derivatives of E^{u) and Ey{u) are typical. 
The next step is to show that the probability density dis- 
tribution P{a) of a at local minima of E vanishes linearly 
with a in the limit a — > 0. Indeed, 

P{a) = {5{E"-a))\^,^, 

^h- j du5{E" ~ a)5{E') ^ \ 
\ ' 
^\amE")5{E')) 

= \amE':,+E'i)5{E'^+E'^)). (79) 

The last equation explicitly demonstrates the aforemen- 
tioned cancellations between the derivatives of E^ and 
i?>. The value of the delta-function product average 



is determined by properties of typical configurations. 
Hence, the low probability of having unusually shallow 
local minima is entirely due to the first factor on the 
left-hand side of Eq. (|79| ) and so -P(a) oc \a\. In partic- 
ular, P c>c e in the case of interest Wm. But this is not 
yet the end of the story. Following Aleiner and RuzirO 
we argue that the probability of having a shallow global 
minimum is additionally suppressed. Indeed, if the co- 
efficient (3 in front of the the cubic term in Eq. (|7^) is 
too large, E{X) would have a second minimum which is 
deeper than that at the ground state X = 0. To avoid 
that \P\ must not exceed ^y/a oc ^/e. The typical depen- 
dence of the resultant E{X) is illustrated in Fig. ||. 



E 




X 

FIG. 8. Energy of the system as a function of the col- 
lective coordinate X. 

The total probability density of the optimal fluctua- 
tion is therefore proportional to e • ^/e ~ with s — 3/2, 
as we claimed. 

It is very likely that s ~ 3/2 is in fact the general re- 
sult independent of the number of dimensions because the 
optimal fluctuation would always have only one soft dir 
rection, i.e., the system is essentially one-dimensionalo 

Finally, let us compare our results with those the 
literature. One group of works, by Feigelman et aZ£j has 
been already mentioned before. We borrowed from them 
the idea of splitting the system into two statistically in- 
dependent parts. The major mistake of these authors 
is overlooking the possibility of frustrations in the sys- 
tem. In other words, they did not realize that it is not 
necessary to require that both £'< and i?> have a small 
second derivative with respect to u when only their sum 
i?< -I- Ey needs to be so. Another point of disagreement 
between us and them is the size of the localized phonons. 
Feigelman et al. take for granted that it is of the order 
of -^//I/e, while we gave an argument that it should be 
Rc, i.e., much smaller. 

In another large group of papers, Refs. p|, [To| , p^ , ^ , 
Vuiiui) OC for the density of states, and a similar depen- 
dence, 'R.euxxi'^) oc w^, for the conductivity were calcu- 
lated. To facilitate the comparison we should point oujt. 
that the conventionally defined phonon density of statesti 
Vuioj) is related to our ^{e) by 

v^{uj) = 2v{puj^)puj. (80) 

Hence, our result is Keaxxii^j) oc ly^iuj) oc cu'^, same as 
in Refs. ^ (see also Ref. ^). All the papers in the sec- 
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ond group are explicitly or implicitly based on the SCBA 
(see Appendix g for more details). The second power of 
Lu can be traced to the unphysical square-root singularity 
v{t) oc |e — ethl^^^ near a band edge eth (cth = in out 
case), which is a well-known basic flaw of the SCBA.E2I 
Note however that even within the SCBA v{€) conforms 
to the generic form (^6|), except it corresponds to s = 1/2 
instead of what we believe is the correct result, s = 3/2. 



VI. QUANTUM AND THERMAL EFFECTS 

So far we have neglected any effects of quantum nature 
or due to a finite temperature. Some of them will be 
addressed in this section but their systematic, in-depth 
treatment is deferred for future work. 

Let us begin with noting that the electrons of the 
pinned WC constantly fluctuate around their equilibrium 
positions. The typical size of such fluctuations is easy to 
find. For example, in strong magnetic fields and at low 
temperatures (fc^T <C hujc) where all the electrons are 
confined to the lowest Landau level, we have 



coth 



1 



2kBT 



n A^(q) + A(q) ^ 
iio H -z q 



n(q) = (pa;,)-V[no +M(q)9')][no + A(q)<z2] 



(81) 
(82) 



This formula can be used away from the thermal or quan- 
tum melting transitions, where the phenomenological cri- 
terion (u^)^/^ < ea is |aatisfied. Here e ^ 0.2 is the 
Lindenmann parameter.EI Incidentally, at T = Hq = 0, 
Eq. (^l]) gives the variance of the electr|©ffl fluctuations in 
the correlated WC of Lam and Girvin.EZi 

Let us focus on the range of temperatures ksT <C htus 
where ujb ^ 'iir fi/rrieUJc is the magnetophonon bandwidth 



(usually, Tiwa/fes ~ 2K). At such temperatures (u 



2)1/2 

is totally dominated by quantum fluctuations and is of 
the order of Ib- We expect that as long as (u^)^/^ is 
smaller than the correlation length of the pinning poten- 
tial ^, which is the smallest length scale in the classical 
theory, all the results obtained in the previous sections 
acquire at most raiaor corrections. j— , 

Recently Fertigia and Chitra et a/.,E3 studied the op- 
posite limit, S, ^ Ib, and found a novel dependence of 
the pinning frequency LUp on the magnetic field, although 
they disagree with each other on the functional form of 
such a dependence. To sort things out we offer a differ- 
ent perspective on this question. Recall that the classical 
theory predicts the l/i3-behavior,Ej 



C(0)c 



(83) 



which follows from Eqs. (y), (g5|), and (|3g|). We pro- 
pose that a reasonably accurate estimate for ujp in the 



regime Ib can be obtained within a quasiclassi- 

cal approximation. The electrons are visualized as com- 
pact wavepackets whose centers of gravity perform clas- 
sical motion, while the quantum effects are presumed 
to be contained in the form- factor of the wavepakets. 
This type of quantum effects amount to replacing the 
bare pinning potential U{r) by its convolution with the 
form- factor Fe{r) of the wavepackets. For small fluctu- 
ations Ib ^ ci the appropriate form- factor is Gaussian, 
Fe(r) = (7r(u^))~-^ exp(— r^/(u^)). Upon the convolu- 
tion, we obtain an effective random potential, with the 
correlation length Is and variance C(0)^^/Z^. In view of 
such modifications, Eq. (p3) transforms into 



(84) 



c(o)c's^ ^<Ib. 



The quasiclassical approximation is certainly not exact; 
nonetheless, when used away from the quantum Hall frac- 
tions, at sufficiently low T, and for robust quantities such 
as LUp, it should be qualitatively correct. 

The linear in _B-dependence of Up derived by Fertig [see 
Eq. (44) in Ref. |l8| is at odds with the i?^-dependence 
found above [Eq. (^4|)] and in Ref. ^ Although Fer- 
tig consider a Poissonian disorder (potential wells or 
"pits" of size s < Ib, energy AV, and the areal den- 
sity Ui < I/ttZ^) rather than the Gaussian one, the real 
source of the discrepancy is his ascertion that the elas- 
tic distortion at the Larkin length scale is of the order 



of n. 



-1/2 



Instead, Ib should be used. Upon this cor- 



rection, the formula for uip acquires the above form ( |84| ) 
with C(0)$2 replaced by n^s'^AV^. 

The finite amplitude of the quantum fluctuations 
has another important repercussion: nonvanishing high- 
order in u terms in the pinning energy. As a result, the 
absorption line acquires extra broadening even at T = 0. 
It originates from finite widths r„ of the energy levels 
n = 1,2,... of the localized magnetophonon oscillators 
(the ground state n = remains unbroadened) . At zero 
T only the transitions between n = and n = 1 levels 
contribute to the absorption. The level width Fi of the 
n = 1 state is determined by the decay of the given mag- 
netophonon into two other magnetophonon excitations of 
lower frequencies, huJi — > hLj2 + ?iw3 (see Fig. |^). 

Such a process is always possible in a truly random 
system, the necessary three-phonon matrix elements, 

M„^^(r) = -neJ2KaKpK^Uir)smK[r-u^''\r)], 

K 

being provided by the cubic anharmonicities. The corre- 
sponding contribution Im H3 to the imaginary part of the 
self-energy can be obtained by evaluating the diagram in 
Fig. |, 



ImH3(q,w) 



2h 
'Z2 



duj' 
27 



k k' 
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X [nsiu}') + n_B(w - uj') + 1] 

X linD^,{k,u}')lmD,r{k' ,uj - uj'), (85) 

^si^) = /, 7- (86) 

At ksT <ti fiLOp only the frequency interval Q < lo' < lo 
contributes to the integral, which leads us to the estimate 



Inin3(ci;p) ^ Imn(ci;p/2 



1% |Imn(wp/2) 



(87) 



This formula is derived assuming that ^ ^ Ib] other- 
wise, in the spirit of the quasiclassical approximation, we 
have to replace ^ by Ib- Even in that case Imn3(a;p) is 
smaller than the previously found Imn(a;p) by a large 
factor ~ {RrJ af'iX/ fiY . At higher temperatures, the 
anharmonic effects are enhanced by a factor of ksT /fnop 
originating from the Bose distribution functions in 
Eq. (^5|). They should start to affect the linewidth when 
kBT ~ {Rc/af^hiop (cf. Ref. H). 



flCL>2 



fid) I 




k', cj -cij 



FIG. 9. Top; low-temperature absorption in a system 
with energy levels broadened by anharmonicities. Bottom: 
the diagram used for estimating Fi. 

A word of caution is in order. The diagram in Fig. ^ 
does not contain vertex corrections, which contain infor- 
mation about, e.g., the soft modes. The soft modes are 
crucial for the dynamical response in strong magnetic 
fields and at the same time are potential sources of much 
larger anharmonicities. Thus, the response of a weakly 
pinned quantum WC may prove to be more nontrivial. 

Concluding this section, we would like to re-emphasize 
that for the problem studied in the main part of the pa- 
per (zero-temperature linear-response of a classical WC), 
the anharmonicities are not important. 



VII. COMPARISON WITH EXPERIMENT 

In this section we will-attempt to analyze the most 
recent experimental datallSO in the light of our under- 



standing of the pinning mode. Let us first discuss the 
position of the pinning resonance. 

From Eqs. ( ^3[ ) and (84) we see that depending on the 
type of disorder, ojp can either decrease as l/jS,or in- 
crease as when the magnetic field increases.ll3 Weak 
magnetic .field dependence of ujp found in the experi- 
mentstSEIl leads us to conclude that ^ remains of the 
order of Ib in the limited range 10-^15 T of the mag- 
netic fields where the pinning mode was detected, i.e., 
^fulsfn 75A. The source of such a short-range random 
potential is likely to be the roughness [©f the heterostruc- 
ture interface as suggested by Fertig.tSI The experimen- 



tal value of 



X 10^ 



for 



5.4 X 10^^ cm- 



can then be reproduced with reasonable values of fitting 
parameters (see Sec. [v| ) s = 30 A, AV = 4.3 K, and 
ni = 4 X lO^^cm"^. Unfortunately, the disorder param- 
eters are poorly known, and so we have to use Up to 
estimate the disorder characteristics, not the other way 
around. For example, we can solve Eq. (Q) for the rms 
amplitude •y/C(0) of the pinning potential, which gives a 
value of the order 0.2 K. 

Let us now discuss the density dependence of ujp. Ex- 
plicitly, it enters Eq. (|8j) only through the shear mod- 
ulus /i(ne)- The latter should behave as cx rie^^ (see 
Sec. O) away from the quantum Hall fractions, e.g., the 
■|-filling. If C(0) is ng-independent, we therefore expect 



iUp (X n^j-i ■ Such a dependence was indeed observed by 
Li et alES for the concentration of holes (p-type samples 
were used) in the range 3-5 x 10^*^ cm^^, which corre- 
sponds to the filling factor range 0.1-0.16. At lower con- 
centrations the pinning frequency continued growing as 
rif, decreased but less rapidly. This comes presumably 
from the fact that the pinning is no longer weak at such 
low Tie', one can verify that Rc is approximately 6a at the 
highest densities, but approaches a (the lattice constant) 
at the lowest densities used in Ref. ^ (In Ref. |lj the 
density dependence of tUp was not investigated) . 

The relationships among the empirical values of ujp 
and Awp and the experimental parameters (rie and B) 
provide another means to test the agreement between 
the theory and the experiment. To this end we rewrite 
Eq. (0) (with s ~ 3/2) in the following way: 



Q = 



AoJr, 



26 



3/2 ' 
eC Tie 

K Bfpk 



1/2 



- 5 



3/2 



(88) 



where fpk = ujp/(2n) is the pinning frequency in cycles 
per second. The largest quality factor Q sa 8 reported by 
Li et al. was achieved for ~ 5.4x 10^" cm^^, _B = 13 T, 
where /p^ was measured to be 1.4 GHz. In Ref. ^ Q « 45 
was found for comparable ne, B, and fpk- Our theoret- 
ical estimate from Eq. (^) is Q sa 260 and far exceeds 
both. We should point out, however, that specific details 
of the experimental setup become important for observ- 
ability of such a high Q. Henceforth we focus on the work 
of Li et aL,E3 since most of their data is available in the 
published record. 
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Because of the finite distance (30 /zm) between the 
plates of the coplanar microwave waveguide (see Fig. 10), 
these were not the ideal q = measurements. We esti- 
mate that the spread of q's imposes un upper bound of 
about 30 for the largest observable Q. The remaining dis- 
crepancy seems to be rooted in finite-temperature and to 
a lesser degree finite incident power effects ignored in the 
derivation of Eq. (|8^). In the experiments, the linewidth 
was decreasing roughly linearly as the temperature varied 
from T = 200 mK to T = 50 mK with other parameters 
held fixed. Below 50 mK the linewidth appeared to reach 
saturation, but only at the nominal level of incident mi- 
crowave power (80 pW). When lower power levels were 
used, the line continued sharpening up. At the lowest ex- 
perimental temperature of 25 mK, the linewidth did not 
show any signs of saturation (this time as a function of 
power) even upon ten-fold input power reduction, which 
was near the sensitivity limit. 



30yum 45jum 



0.6//m 



30/im 



2D WC 





FIG. 10. Top: schematics of the microwave setup used 
by Li et al. The transmission line is shaded, the 2D elec- 
tron (or more precisely, hole) gas is represented by the 
thick horizontal line underneath (adapted from Ref. p^ ). 
Bottom: typical formi of the electric field distribution in 
the plane of the WC.e3 Solid line is the signal in phase with 
the source, dashed-dotted line is the signal in quadrature. 



Several sources of the thermal broadening can be envi- 
sioned: (aV-thermally excited single-particle excitations 
of the WC,I13 (b) phonons of the host semiconductor, (c) 
the anharmonisms of the coUcGtive modes of the pinned 
WC, and perhaps some others .tS 

The classical activation energy of vacancies and inter- 
stitials is of the order of 0.5 K at the densities studied, 
and mechanism (a) should be frozen out at 25 mK. Fur- 
thermore, for the WC of low density the bandwidth of 
such excitations is very narrow and so they are easily 
localized by the random potential. 

The weak electron-phonon coupling in GaAs and small 
phonon phase space at the frequencies involved (ujp) are 
likely to render the mechanism (b) inefficient. 

The preliminary estimate of the broadening due to 
mechanism (c) was given in the previous section. It is 
too small to explain the experimental observations. We 
speculate that stronger anharmonic effects and a better 
agreement with the experiments may be obtained if one 



properly accounts for untypically large anharmonicities 
at the locations of the soft modes, which as we showed 
in this paper, play a very important role in the response. 

It is also quite possible that the actual pinning poten- 
tial is more complicated than the one studied here. In 
reality, it pCfin be due to a combination of the interface 
roughnessKLand dilute residual ions in the vicinity of the 
WC plane.CJ In this case ujp may be determined by the 
former, while Q could be limited by the latter. Such more 
complicated models as well as the dynamic response of 
a pinned quantum WC are interesting subjects awaiting 
further investigation. 

It would be also interesting to investigate if the dis- 
cussed phenomena appear in the conventional charge- 
density waves (see Sec. |). These materials are three- 
dimensional, and some important modifications of the 
present theory may be needed. There are other compli- 
cations such as an extra dissipation due to uncondensed 
quasiparticles but those can be suppressed by lowering 
the temperature. In any case, it would be remarkable if 
narrow pinning modes could be produced in these ma- 
terials simply by applying a sufficiently strong magnetic 
field. 

Our theory should literally apply to the pinning mode 
of a WC formed by electrons on solid hydrogen (see the 
second book cited under Ref. H). We hope that this pa- 
per will encourage further experiments on these or other 
numerous systems where the pinning manifests itself. 
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APPENDIX A: MODEL CALCULATION FOR L-L 
SCATTERING 

In order to evaluate the importance of the L-L scatter- 
ing channel, let us consider the model action 

Me)- I I ul[ci)[Hl{ci) - e]uM) 



ul{v)SL{v)uL{r). 



(Al) 



where Sl is a local operator of a Gaussian white-noise 
type with parameters 

{Sl) - Ho, {Sl{v^)Sl{v2)) = - F"(0)5(ri - v^). 
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Our objective is to calculate the disorder-averaged 
Green's function 



^L(q,e) 



L2 



/ VuLVuluL{(i)ul{(i)e 



where an infinitesimal imaginary correction is assumed 
to be included via e ^ e + iO to make the integrals con- 
vergent. |— I 

Using the "supersymmetry" technique Jill Dl can be 
represented in the form 



'D(j)D(j)*ULU*^e 



where 4>\ — [u^ v^] is a supervector, ?;i(q, w) is an aux- 
iliary fermionic field, and ^^[0] is the sypersymmetric 
Eucledian action 

As = \j 0t(q)(yg + no-e)0(q) 



Consider the effect of such an integration on first. 
Initially, gji remains small and we need to take into ac- 
count only the three lowest order in diagrams (see 
Fig. |l|a) to find 



9r{1 + 51) ^ gnil) + -gl 



Mr 



(27r) 

A<q<A 

t) + (s ^ u) 



Mr, 



(A4) 



where = = (qi - q2)^, T = = (qi - qs)^, and 
U = — (qi — q4)^ are the Mandelstam variables with 
cii, i — I, . . . ,A being the incoming momenta. Similarly, 
the tadpole diagram shown in Fig. |ll[ b determines the 
variation of the renormalized mass Mji, 



Mr{1 + SI) 



1 



' A<q<A 



(27r)2 |q|+M^ 



(A5) 



+ o^"(0) / [0^(r)0(r)] = 



(A2) 



Note that the coefficient in front of the quartic term is 
imaginary, which makes the action non-Hermitian. The 
sign of the imaginary part is such that excitations above 
the vacuum state = decay after travelling a certain 
distance. This distance is simply the phonon mean free 
path. 

It is easy to establish some of the properties of the 
propagator right away. Since the field (f) is massive 
for e < Ho at the tree level, the imaginary part of the 
propagator Dl is small at such e. However, IvuDl should 
grow rapidly as soon as e crosses the e = IIq threshold. 
Below we are going to verify this by explicit calculations. 
Although our methods do not work in the immediate 
vicinity of e = IIq, the expressions obtained for e < IIq 
and e > Ho can be smoothly matched onto each other. 

Rescaling the fields in Eq. (^) by ^J2/Y, we obtain 
the action with the Lagrangian 



C 



M 



(A3) 



where K is the operator, which corresponds to multi- 
plication by |q| in the momentum representation, M = 
2(no-e)/r plays the role of mass, and g = -2V"{0)/Y^ 
is the bare coupling constant. The coupling constant 
is dimensionless, and it is possible to show that the 
field theory is renormalizable. Let us derive the Wilson- 
style renormalization group (RG) equations for the renor- 
malized parameters gn and Mr. As usual, it is done 
by successive integrations over narrow momentum shells 
A{1 + 51) < q < A{1) where A{1) = A^e"' and A^ 1/a 
are the running and the bare high-momentum cutoffs, 
respectively. 



(a) 



FIG. 11. Diagrams for the 
(b) mass term. 



coupling constant and 



In addition to gn and Mji, the kinetic term K also gets 
renormalized. We will neglect such an effect because it is 
of higher (second) order in gji. Consequently, the spec- 
tral function Im Dl is determined just by the mass. 



Im DL{q, e) 



ImMii 



Y {q + ReMi?.)2 + (ImMfl,)2 ' 



(A6) 



In this equation Mji = Mji{q,t) stands for the renor- 
malized mass at such large negative I where the RG 
flow eventually stops either because q ^ A{1) or because 
\Mii\ A (Z). Below we consider exclusively q — case. 

From Eq. ( |A4| ) we find that in the limit jq^l, \Mji\ < A 
function /3 = dga/dl depends only on gji, 



1 



/3 = -9r, 

TT 

and so the RG flow equation is easy to solve: 



9R 



1 - 

9 



(A7) 



(A8) 



At I « Il—t^, 9r becomes of the order of one, and the the- 
ory enters the strong coupling regime. The corresponding 
spatial scale is Ll-l A~^ exp{TT/g). Since the strength 
of coupling maps to the strength of the effective disorder 
in the original formulation, L^-l has the meaning of 
the L-phonon localization length (due to L-L scattering). 
Note that gn seems to exhibit a divergence (Landau pole) 
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at I = II. This is an artifact of the one-loop RG. We ex- 
pect that instead gn flows towards its fixed point value 
of the order of one. 

The ot her scaling equation, for Mn, can be deduced 
from Eq. ( |A5|) , 



7m 



dl An A(0 + Mr 



i5' 



(A9) 



The solution cannot be not expressed in elementary func- 
tions but we can describe its asymptotics. 

1. ((7/47r)A, = < \M\ < A^. 

Here the renormalization of the coupling constant is 
not important because the RG flow is effectively cut off at 
I ^ ln(Ac/|Mji|) before gn manages to get large. Equa- 
tion (A9) reduces to the SCBA equation 



Mr, = M- 



dqq 



An J q + Mr - iS ' 

A 

with the approximate solution 

Re Mr ~ M ~M^, 

IuiMr ~ (5/4)e(-ReMi^)ReMi^, 



(AlO) 



(Alia) 
(AUb) 



where Q{x) is the Heaviside step-function. A word of 
caution is in order here. The imaginary part of Mr, van- 
ishes at positive ReM^ only at the level of the one-loop 
RG. Nonperturbative methods indicate that IbiMr is 
nonzero albeit exponentially small (see below). 



2. Llt^ < \M-M,\^M, 



Now instead of Eq. ( |All[ ) we have 



Im Mr 

Re Mr 
= In ( — 



AIl-L 



-e{-ReMR)ReMR, 



Ir 



Ac 
Re Mr 



1/4 



(M-M,), 



(A12a) 
(A12b) 
(A12c) 



The above comment about small but nonzero ImMj^ 
at ReMii; > applies here as well. 



3. \M-M,\<Lll^. 



Since the one-loop RG cannot be trusted beyong the 
point I = II — TT, we choose to terminate the RG proce- 
dure at this stage. Unfortunately, this leaves us with a 
theory which (a) does not have any small parameters and 
(b) includes not only the original quartic but also higher 
order interaction terms generated by the RG. However, 



on physical grounds we can expect that in this, essen- 
tially massless, limit, the mass scale must be set by the 
running cutoff A{Il), i.e.. 



ImMfl - Acg 



-l/4g-7r/g^ 



(A13) 



This can be c ompared with the solution of the SCBA 
equation (AlO) in the massless limit (ReAf^ = 0), 



ImM|CBA 



A, 



-4ir/g 



(AM) 



Clearly, the SCBA misses the factor of four in the expo- 
nential. 

As mentio ned ab ove, RG-enh anced perturbation the- 
ory formulas ( Allb| ) and ( Al2a) do not work at low ener- 
gies. Nonvanishing Im Mr can be detected usirw-jthe non- 
perturbative method of constrained instantons.ES The de- 
tails of the calculation will be presented elsewhere. Here 
we only quote the result. 



Im Af, 



47re^ 



R 



^ In ^ ^"^^ 



9gl \ReMR 



Ac 



(A15) 



which is va lid for ReM R ^ gR <C 1. Equa- 

tions ( [Al^ ) and ( |AT^ ) match for \ReMR\ - 1/Ll-l- 
Therefore, we expect that the spectral function ImD^ ex- 
hibits a quasi-Lorentzian peak centered at Cp = YM^ ~ 
Hq of exponentially small width 



Aep = I ImA//fl(0, e^) ~ AcYg-'^e~^^3 _ 



(A16) 



APPENDIX B: DYNAMICAL RESPONSE 
WITHIN THE SCBA 

Within the SCBA the self-energy Hapiuj, q) is diagonal, 
Ilap(oJ,q) = SapTl and weakly g-dependent for qa <C 1. 
H is to be found from the self-consistency equation ( |37| ) 
which we reproduce here for convenience, 



V"{0) J tTD{k,uj). 



(Bl) 



As discussed in Sees. ^ and the SCBA applies only 
at relatively high frequencies. Moreover, in strong mag- 
netic fields and near the pinning frequency lo ujp, it 
is not even qualitatively correct. One may ask why we 
bother elaborating on this faulty approximation scheme 
here. The reason is as follows. We discovered that as 
far as the dynamical response .is .pnncjerned, all alterna- 
tive theories suggested so farESEjIla'E^ are merely differ- 
ent forms of the SCBA. Thus, we deemed that it would 
be helpful to expose their interrelationship and correct a 
few calculational errors. 

Let us briefly summarize the resul1js-|0btained by pre- 
vious authors. Fukuyama and LeeEiJ concluded that 
Aojp ~ ujp but did not present the details of the cal- 
culation. Perturbative analysis of FertigpS yields an ex- 
ponentially small AiVp/up (for weak disorder). Gaussian 
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variational replica method (GVM) of Chitra et al^ re- 
duces to the SCBA at nonzero frequencies, and when an- 
alyzed further, predicts a power law dependence of Aujp 
on the disorder strength. Finally, the SCBA analysis 
in the Appendix of Ref. ^ leads to yet another depen- 
dence. One reason for such a disarray of conflicting re- 
sults is the extreme sensitivity of the solution 11 (cj) of 
Eq. (Bl) to the precise value of the dimensionless param- 
eter C = -F"(0)/(47r/ino), where Uq = U{0). In prin- 
ciple, C is fully determined by 5*0, but 5*0 is known only 
approximately. A strong partial cancellation of 5*0 by the 
second term in Eq. (Bl ) leaves us only the order of magni- 
tude estimate C ~ 1. One could argue, as Fukuyama and 
Lee did in analogous situation,u that C = 1 is the "best" 
choice because others lead to various physical absurdities. 
For example, if C > 1, then Tl{uj) acquires a finite imag- 
inary part at complex to, which means that the ground 
state is unstable; if C < 1, then n(w) has an imaginary 
part only above some threshold frequency Wth > 0, i.e., 
the phonon spectrum has a gap, which also appears to 
be unphysical. The truth is, of course, that the SCBA 
is simply unable to correctly describe the properties of 
the Lo tail; therefore, it cannot be relied upon for 
selecting a "good" value of C. More sensible approach is 
to investigate a certain range of C around unity, hoping 
that certain quantities depend on C only weakly. Alas, 
the linewidth Aujp is not one of such quantities. 

We found that the results of Chitra et al. are recovered 
for the "Fukuyama-Lee choice," C = 1. (Therefore, these 
two groups of authors should have obtained the same re- 
sults). A slight reduction of C from unity is sufficient to 
cross over to very different predictions of Fertig. Let us 
now demonstrate this in more detail. 



Using the definition of C and Eq. (Bl), we obtain 



n(cj) - Ho = 2^Cno / dkk [Drik, 0) - Drik, oj)] . (B2) 



Let us now describe the solution of this equation 

for different C in the weak pinning regime Rc S> a, which 
corresponds to the inequality a = ^/TToJl/Y = /i/A <C 1. 



leads to Awp/wp a 



1/2 = (^/A)i/2. Although this re- 
sult is of the same form as Eq. (|66|), with s = 1/2, we 
believe that s = 3/2 is correct, see Sees. IV and M. 



2. 1 - C> a. 



Equation (B3) still applies as long as LUpQ — u is positive 
and not too small. At larger frequencies, 11 (w) receives 
additional contribution from the pole fc, of Dt in the 
complex plane of fc, which moves close to the real axis. 
The resonance linewidth Aujp is exponentially small, 



Imn(u;p) 



pUJc 



(B4) 



where F ~ n(tJth)/(Cno) — 1. The threshold frequency 
LUth is only slightly below Wp, by an amount of the order 
of Aujp. To be exact Wth can be found from the condition 



that the solution of Eq. (B2) is also the solution of the 
"derivative" of this equation with respect to 11, 



1 2/xCno 



dkk[eec + (n + Yk)^] 
[{U + fik^){U + Yk)- eec 



(B5) 



At frequencies somewhat higher than LUp, Imn(w) is dom- 
inated by the aforementioned pole at fc, = {eec — lP)/UY 
near the real axis, |Re fc* | ^^Im fc, | . Im 11 is given simply 
by the residue of this polejla 



Imn(cj) ~ 27ra^Cno 1^-1 



(B6) 



In conclusion, we would like to reiterate that neither of 



Eqs. (|B3D , (B4), or (B6) describes the dynamical re- 
sponse of the system correctly. We derived them here 
just to facilitate the comparison with the previous work 
on the subject. The correct expressions for Imll and axx 
are given in Sec. IV. 



1. < 1 - C < a. 
In this case 

n(c^) ~ CTo + mo ^{uj/[bp)^~ii~c)^, (B3) 

where uip = ujp/{TTay^^ and ujp ~ Ylo/pujc- As one can 
see, the threshold frequency is wth = (1 — C)ojp. For 
C = I Eq. (p3|) coincides with that of Ref. |l^. Strictly 
speaking, Eq. (|B3|) describes the solution of Eq. (^ ) only 
for uj ujp. Near the resonance, lu cUp, it is off by a 
numerical factor of the order of unity. In principle, this 
factor can also be calculated. We limit ourselves to the 
order of magnitude estimate Imn(wp) ^ a^/^IIo, which 
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